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Abstract. We present empirical fits to the UBVRI light 
curves of type la supernovae. These fits are used to objectively 
evaluate light curve parameters. We find that the relative times 
of maximum light in the filter passbands are very similar for 
most objects. Surprisingly the maximum at longer wavelengths 
is reached earlier than in the B and V light curves. This clearly 
demonstrates the complicated nature of the supernova emis- 
sion. 

Bolometric light curves for a small sample of well- 
observed SNe la are constructed by integration over the opti- 
cal filters. In most objects a plateau or inflection is observed in 
the light curve about 20—40 days after bolometric maximum. 
The strength of this plateau varies considerably among the in- 
dividual objects in the sample. Furthermore the rise times show 
a range of several days for the few objects which have obser- 
vations early enough for such an analysis. On the other hand, 
the decline rate between 50 and 80 days past maximum is re- 
markably similar for all objects, with the notable exception of 
SN 1991bg. The similar late decline rates for the supernovae 
indicate that the energy release at late times is very uniform; 
the differences at early times is likely due to the radiation dif- 
fusing out of the ejecta. 

With the exception of SN 1991bg, the range of absolute 
bolometric luminosities of SNe la is found to be at least a factor 
of 2.5. The nickel masses derived from this estimate range from 
0.4 to 1.1 Mq. It seems impossible to explain such a mass 
range by a single explosion mechanism, especially since the 
rate of 7— ray escape at late phases seems to be very uniform. 

Key words: supernovae: general - Stars: fundamental parame- 
ters 



1. Introduction 

The temporal evolution of a supernova's luminosity contains 
important information on the physical processes driving the ex- 
plosion. The peak luminosity of a Type la Supernova (SN la) 
is directly linked to the amount of radioactive 56 Ni produced 



Hoflich et al. 1997, Eastman 1997, Pinto & Eastman 2000a) 



in the explosion (Arnett et al. 1985, Branch & Tammann 1992, 



The rise time of the light curve is determined primarily by 
the explosion energy and the manner in which the ejecta be- 
come optically thin to thermalized radiation, i.e. the opacity 



(Khokhlov et al. 1993). The late decline of the light curve is 



governed by the combination of the energy input by the ra- 
dioactive material and the rate at which this input energy is 
conve rted to optical photons in the ejecta (Leibundgut & Pinto 



1992) 



The apparent uniformity of SN la light curves in photo- 



graphic (pg), B, and V filters (Minkowski 1964) prompted 
the adop tion of standard lig ht curv e templates (e.g. Elias et 
al 



1985 



1991b 



Dogg ett & Branch 1985 , Leibundgut et al 
Schlegel 1995 ). Early indications that the standard templates 
fail to describe the full range of SN la light curves came from 
the observations of SN 1986G, which displayed a much more 
rapid evolutio n than any other SN observed up to that time 
(Phillips et al. 1987 1. The demise of the simple standard candle 
treatment was brought about by the observations of the faint 



SN 1991bg (Filippenko et al. 1992, Leibundgut et al. 1993 



Tu- 



ratto et al. 1996[ ) and the subsequent derivation of a correla 



tion between peak luminosity and decline rate after maximum 
(Phillips |1993| , Hamuy et al. |l996a| , Riess et al. |l996a| ). A clear 
demonstration that SNe la do not display a uniform photomet- 
ric evolution was provided earlier from the infrared J, H, and 
K light curves (Elias et al. 1985, Frogel et al. 1987). Observa- 



tions and analyses of the near-IR R and / light curves (Suntzeff 



1996, Vacca & Leibundgut 1997) confirm this variation. The 



red and near-infrared light curves exhibit a second maximum ~ 
20 to 30 days after the B peak. This second maximum occurs 
at different phases and with differing strengths in individual 
SNe la and in at least one case (SN 1991bg; Filippenko et al. 



1992, Leibundgut et al. 1993, Turatto et al. 1996) is altogether 



absent. 

The light curves of SNe la are often described as a one- 
parameter family. The correlation of the BVRI light curve 
shapes with the peak absolute magnitudes has been employed 
to improve the distance measurements derived from SNe la 
(Hamuy et al. |1996a| , Riess et al |1996aj Garnavich et al. |1998 , 
Schmidt et al. I1998L Riess et al. |1998a|) and is of fundamental 
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importance to keep the systematic uncertainties in the deriva- 
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tion of cosmological parameters small. Techniques that fit stan- 
dard templates (e.g., Hamuy et al. 1996a), or modified versions 



of these templates (e.g., Riess et al. 1996a; Perlmutter et al 



1997), to the observed light curves make use of this one pa- 



rameter description of the light curves. These methods have the 
advantage that they can be applied even to rather sparsely sam- 
pled light curves. In a more extreme form, the photometry can 
be supplemented by spectroscopy to provide a distance mea- 
surement with minimal data coverage . This " snapshot" method 
has been advocated by Riess et al. ( 1998b ). All these meth- 
ods make the assumption that SNe la form an "ordered class". 
This description is validated by the improvement in the scatter 
around the linear expansion line in the local Universe and also 
by the fact that new objects can be successfully corrected with 
the correlations derived from an independent sample. 

While clearly useful for comparing local SNe with high 
redshift SNe, the template method does not allow one to in- 
vestigate finer and more individual features in a large sample 
of SN la light curves. Hence, the detailed study of the explo- 
sion and radiation physics cannot be carried out with such an 
analysis. For data sets which are densely sampled, however, 
template methods are not necessary. Recent bright SNe la have 
been observed extensively and very detailed, and accurate light 
curves have become available. Most of these supernovae have 
been used as the defining objects for the templates to correct 
other, more sparsely observed, SNe la. 

To analyze light curves of many SNe la in an individual 
fashion a parameter fitting method which can be applied to sin- 
gle filter light curves has been devised (Vacca & Leibundgut 
1996, 1997). The photometric data are approximated by a 



smooth fitting function. We are not fitting model light curves 
based on explosion physics, but simply attempt to match the 
data in an objective way. We have investigated well-observed 
SNe in a small sample to check our method. It allows us to 
search for correlations among various light curve parameters, 
to accurately fit the filter light curves at any phase, and to con- 
struct continuous bolometric light curves. The application to 
the light curves of SN 1994D has already provided one of the 
first b olometric light curves of SNe la (Vacca & Leibundgut 



1996) 



Observationally-derived bolometric light curves provide a 
measure of the total output of converted radiation of SNe la, 
and therefore serve as a crucial link to theoretical models and 
calculations of SN explosions and evolution. The total lumi- 
nosity from a SN la is much easier to calculate from theoretical 
models than the individual filter light curves, which are dom- 
inated by line blending effects requiring complicated multi- 



group calculations (Leibundgut & Pinto 1992, Eastman 1997 



Hoflich et al. 1997). In addition, the observed bolometric peak 



luminosity of SNe la provides a measure of the total amount 
of nickel synthesized in the explosion and can be used to test 
various explosion models. Although no SN la has ever been 
observed in every region of the electromagnetic spectrum si- 
multaneously, fortunately bolometric light curves can becon- 
structed al most en tirely from optical data alone (Suntzeff 1996 , 
Leibundgut |l996 l. 



In Sect, g we describe the sample of objects and the data 
set used. The parameter fitting method employed for our anal- 
ysis is described in Sect. |3|. This is followed by a discussion of 
the phases of the maximum epoch in different filters (Sect. Q). 
The construction of the bolometric light curves is presented in 
Sect, ^followed by a discussion of the uncertainties in the bolo- 
metric light curves. We summarize our findings in Sect. ^ our 
conclusions are given in Sect. [7| 

2. Observational Data 

The data for our analysis come from recent photometric obser- 
vations of bright SNe la, consisting of the large data collections 
presented by Hamuy et al. ( |1996b ) and Riess (1996 see also 
Riess et al. |1999[ ). Table [l] summarizes the objects which have 
sufficient data, in at least the B VRI filters, to construct accu- 
rate individual filter light curves. Listed are the available filters 
(column 2) and the references of the SN photometry (column 
3), the adopted distance modulus (column 4) and its source 
(column 5). The distance modulus is taken either from direct 
distance measurements via cepheids, surface brightness fluc- 
tuations or Tully-Fisher luminosity-line width relation or by 
adopting a Hq of 65 km s _1 Mpc -1 for galaxies in the Hubble 
flow {v > 2500 km s _1 ). We use t he host galaxy dust extinc- 
tion values of Phillips et al. ( 1999 ) (column 6). The Galactic 
extinction given by the COBE dust maps is provided in column 
7 (Schlegel et al. |1998| ). 

Not all objects have sufficient observations in all filters to 
guarantee that our method will work. In some cases, the light 
curve fits had to be restricted to a limited phase range. Only 
three objects (SN 1989B, SN 1991T, and SN 1994D) have a 
significant number of U filter observations. Our analysis does 
not include any JHK or ultraviolet data. 

3. Fitting filter light curves 

3.1. Method 

The light curves are analyzed using a descriptive model (Vacca 
& Leibundgut 1996). For each supernova the observed light 
curve in each filter is fit with an empirical model consisting of 
a Gaussian (for the peak phase) atop a linear decay (late-time 
decline), a second Gaussian (to model the secondary maximum 
in the V, R, and / band light curves), and an exponentially 
rising function (for the pre-maximum segment). The functional 
form of the fit is: 



/o+7(i-<o)+3oexp(^2)!) 



+ gi exp 



-(t-tiy 

2a? 



l-exp(^) 

The first Gaussian and the decline are normalized to the 
phase to, while the second Gaussian occurs at a later phase t±. 
The exponential cutoff function for the rise has a characteristic 
time 6 and a separate phase zero-point r. The amplitudes of the 
two Gaussians, go and g\, as well as the intercept of the line, 
/o, the slope, 7, the phases to an d t-y, and the widths, <jq and 
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Table 1. Photometry of well-observed SN la 



SN 



(1) 



Filter 



(2) 



ref." 



(3) 



DM 



(4) 



ref. a E(B - V) 



(5) 



Phillips et al. ( |1999j ) 

(6) 



Schlegel et al. ( |1998| ) 

(8) 



SN1989B 

SN1991T 

SN1991bg 

SN1992A 

SN1992bc 

SN1992bo 

SN1994D 

SN1994ae 

SN1995D 



UBVRI 
UBVRI 
BVRI 
BVRI 
BVRI 
BVRI 
UBVRI 
BVRI 



4,5 



1 

2 
3, 
6 
7 
7 

8, 9, 10, 11 
12 



BVRI 12, 13 



30.22 
31.07 
31.26 
31.34 
34.82 
34.63 
30.68 
31.86 
32.71 



14 
15 
16 
17 
7 
7 
9 
18 
12 



0.340 
0.140 
0.030 
0.000 
0.000 
0.000 
0.000 
0.120 
0.040 



0.032 
0.022 
0.040 
0.017 
0.022 
0.027 
0.022 
0.031 
0.058 



"References: 
1 - Wells et al 
4 - Leibundgut et al 



1994,2 



1993 



Lira et al 
5 



1998, 3 



Filip penko et al 
Turatto et al. |1996| 6 - Suntzeff 
Richmond et al. 119951, 9 - Patat et al 



1992 



1996 



7 - Hamuy et al. |1996H 8 - Richmond et al. |1995| , 9 - Patat et al. |1996 
10 - Meikle et al. |1996j 1 1 - Smith et al. p riv. co mm., 12 - Riess et al. 
13 - Sadakane et al. |1996j 14 - Saha et al. |l999| , 15 - Fisher e 
16 - Hamuy et al. |l996a| and references therein, 17 - Suntzeff et al 
18 - Riess et al. 1996a and references therein 



1999 



1999 



1999 and references therein, 



<7i, the characteristic rise time 9 and its phase r are free pa- 
rameters in the fit. Each filter light curve is fitted individually 
and independently using a ^-minimization procedure to de- 
termine the best-fit values of the parameters. Other parameters 
which characterize the shape of the light curve, such as the time 
of maximum brightness or Amis, can then be derived from 
the best-fit model. Although the model is a completely empiri- 
cal description of the general shape of SNe la light curves, we 
note that theoretical models (Pinto & Eastman 2000a) predict 
a Gaussian shape for the peak in models with constant opacity 
and 56 Ni buried well within the ejecta. 

Fig. |] shows the example of the fit to the R light curve of 
SN 1992bc. The various components of the fit are displayed. 
The exponential function (solid line) rises steeply to unity, 
modeling the rapid brightness increase of the supernova. The 
decline rate (dotted line) is set by the long tail beginning about 
50 days past maximum. The first Gaussian (dashed line) fully 
describes the maximum phase, as further demonstrated in the 
inset. The second Gaussian (dashed-dotted line) reproduces the 
"bump" in the light curve. The bottom panel of Fig. [I] displays 
the residuals of the fit with the observational error bars. It is 
clear that the function is an accurate, continuous description 
of the data. The small systematic undulations of the residu- 
als indicate that the fit is not perfect and can be used to make 
detailed comparisons between individual supernovae. We have 
fitted this model to the filter light curves of more than 50 su- 
pernovae. These fits will be presented in forthcoming papers. 

The fits produce objective measures of the magnitude and 
date of maximum, the extent of the peak phase, the amplitude 
and extent of the secondary peak, the late decline rate, and an 
estimate of the rise time. The accuracy of these parameters 
depends strongly on the number and quality of the observa- 
tions in each phase. While in most cases the late decline can 
be derived fairly easily, the rise time is undetermined when no 



/ Am 15 






1 5 days : 




1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



8940 8960 8980 
t (JD 2440000 + ) 



Fig. 1. Fit of Eq. (|lj) to the R photometry of SN 1992bc. The 
various components of the fit are displayed as the following: a 
linear (in magnitudes) decay (dotted line), an exponential rise 
factor (solid line) from to 1 as indicated by the right ordi- 
nate, two Gaussians for the peak and the second maximum 
(dashed and dashed-dotted line) with the amplitude in magni- 
tudes as indicated by the right ordinate. The inset shows the 
light curve around maximum light. Amis can be obtained eas- 
ily because of the continuous representation of the data by the 
fit. The residuals from the fit are displayed in the bottom panel 
with the observational error bars. 



pre-maximum observations are available. It is possible to de- 
rive other light curve shape descriptions, such as Amis, for 
each filter objectively and independently without comparison 
to templates. 
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Fig. 3. Distributions of the differences between the best fit and 
the Monte Carlo simulations for the time of B maximum (left) 
9000 9020 9()10 and B magnitude at maximum (right) for SN 1992bc. 



Fig. 2. Fit of Eq. (|IJ) to the observed / data of SN 1992bc. The 
description of the lines is given in Fig. [l] 



Other functional forms could be imagined for the fit. How- 
ever, the late decline (linear in magnitudes) and the steep rise 
can be assumed to have simple functional forms. Fitting poly- 
nomials or spline functions, for example, could not produce the 
transition to the linear decline as observed, but would simply 
transfer the undulations of the early phases to the late decline 
and fail to match the observed linear decline. Furthermore, such 
functions would not provide a small set of adjustable parame- 
ters which can be compared between different objects. 

Nevertheless, the adopted functional form is not a perfect 
representation of all optical filter light curves. When fitting the 
/ light curves of some SNe la we found that the model had 
difficulties to match the observations with the same accuracy 
as in the other filters. This can be seen in Fig. [| The slope of 
the / band light curve decline (dotted line) is so steep that the 
model cannot reproduce the observed light curve. This steep 
slope causes the function to overshoot the observed values near 
the second maximum and the fitting program inverts the second 
Gaussian to produce the observed minimum in the light curve. 
If the fitting routine is constrained to fit only positive Gaus- 
sians, the exponential rising function (solid line) is shifted to 
late times and depresses the linear decay (dotted line) before 
the second maximum; in addition the first Gaussian (dashed 
line) is enhanced and shifted to earlier times. This implies that, 
for the / band, the fitted parameters describing the first Gaus- 
sian (i.e. its magnitude, its center and its standard deviation) 
cannot be used to characterize the first maximum, even if the 
light curve appears fitted well. In spite of this problem, the I 
band light curves of the SNe in our sample can all be reason- 
ably well fit by the model and the data can be represented by 
the continuous fits. This also means that the derived light curve 
parameters, e.g. decline rates, are meaningful. 

This steep decline in the / band is a rather unexpected re- 
sult. The physical explanation behind this behavior is that the 
/ filter light curve is dominated by a rapidly decreasing flux 



component and significant flux redistribution takes place in the 
evolution. Similar results have been found by Suntzeff ( 1996 ) 
in the observations of SN 1992 A and Pinto & Eastman (|2000b, 



see also Eastman 1997 1 in theoretical models. 



3.2. Uncertainties in the model parameters 

The fitting procedure provides an estimate of the goodness-of- 
fit, as well as the associated uncertainties on the fit parame- 
ters. Uncertainties on the derived quantities depend strongly 
on the quality and number of data points. If there are no data 
before 5 days prior to maximum, e.g., the rise time cannot be 
determined reliably. We have estimated the uncertainties us- 
ing a Monte Carlo method. We constructed synthetic data sets 
with the same temporal sampling as the observed light curve. 
The magnitude at each point was computed from the best fit 
model, with the assumption of a Gaussian probability distribu- 
tion whose width was given by the observational uncertainty. 
The standard deviations of the observed data points range from 

0. 03 to 0.14 magnitudes. In this manner, 2000 synthetic data 
sets were simulated. Each synthetic data set was fit and the fre- 
quency distribution for each model parameter was constructed. 
Two examples of the resulting distributions (the time of B max- 
imum and the B magnitude at maximum for SN 1992bc) are 
shown as histograms in Fig. ||. The differences between the 
mean and the best fit is less than 0.03 days for the time of max- 
imum and less than 0.01 magnitudes for the peak brightness. 
In this example the skewness and kurtosis are negligible. Non- 
zero values of skewness and kurtosis provide an indication of 
the unreliability of the derived fit parameters. 

Occasionally the standard deviations obtained for a fit pa- 
rameter from the fitting procedure differed significantly from 
the values obtained from the Monte Carlo simulations. This 
arises because the standard deviations derived from the fit are 
calculated from the diagonal elements of the curvature matrix, 
while the standard deviations derived from the Monte Carlo 
simulations result from actual frequency distributions. (Both 
standard deviations should be the same for a perfect light curve, 

1. e. well sampled data, where the brightness differs from the 
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Fig. 4. Histograms showing the times of the filter maxima from 
U to I and of the bolometric luminosity relative to the B max- 
imum. The vertical lines show the time differences as expected 
for an expanding, adiabatic cooling sphere according to Arnett 
fll982|). 



analytical shape of the light curve by only a small error. We 
confirmed our code with this test.) In all cases we adopted the 
uncertainties given by the Monte Carlo simulations. 



4. Epochs of maxima 

For all 22 SNe la whose light curves were well-sampled 
around the peak (SNe 1989B, 1991T, 1991bg, 1992A, 1992bc, 
1992bo, 1994D, 1995D, see Table 0; SN 1990N, Lira et al. 
19981 ; SNe 1990af, 1992P, 1992al, 1992 bh, 199 2bp, 1992bs, 
1993H, 1993Q, 1993ag, Hamuy et al. |l996b|; SNe 1994S, 



1994ae, 1995E, 1995ac, 1995al, 1995bd, 1996X, Riess 1996, 



Riess et al. 1999) we determined the epochs of individual fil- 
ter maxima. A general trend can be observed in Fig. [I| which 
presents the time of the filter maxima relative to the epoch of 
the maximum in B. The U band appears to rise faster than the 
B band, although the limited number of objects (4 supernovae) 
precludes any definitive statements. The V light curve clearly 
peaks some time after the B maximum for all observed SNe la. 
This behavior is more or less expected from models of an ex- 



panding cooling atmosphere (Arnett 1982), as indicated by the 
lines in Fig. Q However, the behavior of the R and I band max- 
ima do not agree with this simple model; for these bands, the 
light curve maxima are reached earlier than in a thermal model. 
The I histogram is very broad but for most objects the rise in I 
is clearly faster than in B. This trend is further continued in the 
JHK light curves which all pea k before B (Elias et al. 1985 , 
Leibundgut 1988| , Meikle 2000 ). The early appearance of the 



peak in the infrared rules out the idea of an expanding, cooling 
sphere. 



5. Bolometric light curves of SNe la 

The flux emitted by a SN la in the UV, optical, and IR wave- 
lengths, the so-called "uvoir bolometric flux", traces the radi- 
ation converted from the radioactive decays of newly synthe- 
sized isotopes. As nearly 80% of the bolometric luminosity of 
a typical SN la is emitted in the range from 3000 to 10000 A 
(Suntzeff |l996| ), the integrated flux in the UBVRI passbands 
provides a reliable measure of the bolometric luminosity and 
therefore represents a physically meaningful quantity. This lu- 
minosity depends directly on the amount of nickel produced in 
the explosion, the energy deposition, and the 7-ray escape, but 
not on the exact wavelengths of the emitted photons. 

We used the fits of the filter light curves in the UBVRI 
passbands to construct an optical bolometric light curve for the 
SNe in our sample. All objects with well-sampled (U)BVRI 
light curves and sufficient coverage from pre-maximum to late 
decline phases have been included. The objects are listed in 
Table [|. To calculate the absolute bolometric luminosities, one 
has to account for reddening and distance moduli; the values 
we adopted were taken from the literature as listed in Table [fl. 
A galactic extincti on law has been employed, as justified by 
Riess et al. ( |l996b| ). 

The bolometric light curves of our sample of SNe la are 
shown in Fig. |[ Only the time range with all available filter 
photometry is plotted. The peak luminosities are clearly dif- 
ferent for the objects in the sample, and these differences are 
larger than the uncertainties in the derivation. The most strik- 
ing feature, however, is the varying strength of the secondary 
shoulder, which stems from the R and I light curves (see also 



Suntzeff 1996) 



The distribution of the epochs of the bolometric peak lu- 
minosity relative to the B maximum is shown at the bottom 
of Fig. U The time of bolometric maximum suffers from the 
additional uncertainty that the contributions from the different 
wavelength regimes change rapidly before and during the peak 
phase. Since we are not including any UV flux in our calcula- 
tions and SNe la become optically thin in the UV around this 
phase, the epoch of the maximum could actually shift to ear- 
lier times than what we measure. This might explain the dis- 
crepancy with the earlier determination of the bolometric max- 
imum for SN 1990N and SN 1992A by Leibundgut ( |l996| ) and 
Suntzeff (1996), respectively, who included IUE and HST mea- 
surements. 

5.7. Uncertainties introduced by the data 
5.1.1. Missing passbands 

Flux outside UBVRI: In our analysis we have neglected 
any flux outside the optical wavelength range. In particular, 
contributions to the bolometric flux from the ultraviolet (be- 
low 3200 A) and the infrared above 1 /im {JHK) regimes 
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Fig. 5. Bolometric luminosities for 9 well observed SNe. Only 
the time range with the maximum number of filters for each su- 
pernova is displayed. The reddening and distance moduli have 
been taken from the literature as listed in Table |l] and correc- 
tions for missing U band are applied as described in Sect. |5.1.l| . 
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Fig. 6. Correction factors for missing passbands. The correc- 
tion factors are obtained by comparing the flux in the passbands 
with the total UB VRI -flux. 



should be considered in the calculations of the bolometric flux. 
Using HST and IUE spec trophotometry for SN 1990N and 
SN 1992 A, Suntzeff (11996b estimated the fraction of bolomet- 



ric luminosity emitted in the UV. He found that the bolometric 
light curve is dominated by the optical flux; the flux in the UV 
below the optical window drops well below 10% before max- 
imum. The JHK evolution was assumed to be similar to that 



presented by Elias et al. (1985); these passbands add at most 
10% at early times and no more than about 15% 80 days after 
maximum. For example, examination of the data of SNe 1980N 
and 198 ID (infrared data from Elias et al. (1981 1 and optical 
from Hamuy et al. (1991 1) shows that not more than 6% of 
the total flux is emitted beyond / until 50 days past maximum, 
when the IR data stop. 

Corrections for passbands missing in the optical range: For 
missing passbands between U and / one can infer correc- 
tions derived from those SNe la which have observations in 
all filters. Fig. || shows the correction factors obtained from 
SN 1994D. A cautionary note is appropriate here: SN 1994D 
displayed some unusual features, in particular a very blue color 
at maximum. 

To estimate the flux corrections we divided the flux in each 
passband by the bolometric flux. Since the filter transmission 
curves do not continuously cover the spectrum (i. e., there are 
gaps between U and B as well as between B and V and over- 
laps between V and R, and R and /) the coaddition suffers 
from the interpolation between these passbands. 

An interesting result from Fig. ^|is the nearly constant bolo- 
metric correction for the V filter. This filter has been used in 
the past a s a sur rogate for bolometric light curves (e.g. Cappel- 
laro et al. 1997 1. We confirm the validity of this assumption for 
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Fig. 7. Residuals for substituting individual passband observa- 
tions of SN 1989B by the corrections derived from SN 1994D. 
The errors never exceed 10% even when more than one filter is 
missing. 



phases between 30 and 110 days after B maximum where the 
overall variation is less than 3%. 

In order to test our procedure, we calculated bolometric 
light curves for the three SNe la which have the full wavelength 
coverage after purposely omitting one or more passbands and 
applying the correction factors we derived from SN 1994D. As 
Fig. demonstrates, the error is less than 10% at all times even 
if more than one filter is missing, although the errors vary con- 
siderably during the peak and the secondary shoulder phases. 
The results of this exercise gave us confidence that we could 
correct the bolometric light curves of the remaining six SNe 
for the missing U band without incurring large errors. 
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5.1.2. Effects of systematic differences in photometry sets on 
bolometric light curves 

Filter light curves from different observatories often show sys- 
tematic differences of a few hundredths of a magnitude. We 
examined the effect such systematic errors might have on our 
bolometric light curves by artificially shifting individual filter 
light curves by 0. 1 magnitudes and recomputing the bolometric 
light curves. In all cases the effect on the bolometric light curve 
is far less than the uncertainties in the distances and the extinc- 
tion of the supernovae. Shifting the R or / light curve by 0. 1 
magnitude we found that the bolometric luminosity changed 
by 2% at maximum and 5% 25 days later (approximately at 
second maximum). For typical systematic uncertainties of 0.03 
magnitudes in all filters, a maximal error of 2 to 3% is incurred. 

We also constructed the bolometric light curve of 
SN 1994D for the individual data sets available (see references 
to Table [l]). The difference never exceeds 4% out to 70 days 
after maximum even though the photometry in the individual 
filters differs up to 0.2 magnitudes in U, 0.1 mag in BRI and 
0.05 mag in V. 

5.2. Uncertainties introduced by external parameters 

While uncertainties in the distance moduli will affect the ab- 
solute luminosities in each passband (as well as the bolomet- 
ric luminosity), the shape of the light curve is unaffected. As 
all distances used here are scaled to a Hubble constant of 
Ho = 65 kins" 1 Mpc -1 , the luminosity differences are af- 
fected only by errors in the determination of the relative dis- 
tance modulus. 

Reddening, however, affects both the absolute luminosity 
and the light curve shape. The influence of extinction changes 
as a function of phase with the changing color of the supernova. 
A color excess of E(B — V) = 0.05 decreases the observed 
bolometric luminosity at t = < m ax(bol) by 15%, while near 
the time of the second maximum in the R and I light curves 
(t = 20 days) the observed bolometric luminosity is reduced 
by 12%. A stronger extinction of E(B — V) = 0.35 reduces 
the observed bolometric luminosity by 67% (56%) at t = i max 
(f = 20 days). 

The uncertainty in the reddening estimate introduces subtle 
additional effects. If 8E(B — V)= 0.02 mag at low redden- 
ings (E(B — V)< 0.05), an additional uncertainty of 5% in 
the bolometric luminosity is introduced. At higher reddening 
values (E(B - V) = 0.3), uncertainties of 6E{B - 10=0.05 
and 0.10 produce changes of 15% and 31% in the bolometric 
luminosity, respectively. 

The decline rate of the bolometric light curve decreases 
with increasing reddening due to the color evolution of 
the supernova and the selective absorption. For SN 1994D 
Ami5(Bol) would evolve linearly from 1.13 at hypothetical 
E(B - V) = to 0.99 at E(B - V) = 0.35. The linearity 
breaks down at about E(B — V) w 0.5. 



-20 20 40 60 80 100 

t (days after maximum) 

Fig. 8. Constructing the bolometric light curves before and af- 
ter fitting the filter light curves of SN 1992bc. The stars indi- 
cate the bolometric luminosity as calculated directly from the 
observations. The correc tion for missing U band is applied as 
described in Sect. |5.1.lj only the epochs with the maximum 
number of filters are shown. The dotted line is the bolometric 
luminosity calculated from the fits to the individual passbands. 



1982 ). This implies that the color evolution influences the 
shape of the filter light curves. A color difference of B — V= 
0.7 for SNe la in the first 15 days results in an increase of 
Ami5(B) by 0.2 x E(B — V) simply due to the color depen- 
dence of the reddening. The increase of Amis for blue filters 
is larger than for the redder passbands. 



5.3. Uncertainties introduced by the method 

The effect of fitting the light curves before constructing the 
bolometric light curves can be seen in Fig. ^[ The bolometric 
flux for SN 1992bc determined in this way is compared against 
the straight integration over the wavelengths of the filter obser- 
vations. The agreement between the two approaches is excel- 
lent and no differences can be observed. This also applies to 
the correction for the missing U filter observations. Deviations 
can be found only for those parameters derived from the light 
curve, which extrapolate far from the observed epochs, e.g. the 
rise time. 

The integration over wavelengths to calculate the bolomet- 
ric light curve can be performed in different ways. A straight 
integration of the broad-band fluxes has to take into account 
the transmission of and any overlaps and gaps between fil- 
ters. We experimented with several integration techniques, but 
found that for all interpolation methods, we reproduced the 
bolometric flux to within 2% at any epoch we considered. This 
has already been shown for SN 1987 A by Suntzeff & Bouchet 



(1990). We chose to calculate the bolometric flux by summing 



As described by Leibundgut ( |1988[ ) the reddening depends 
on the intrinsic color of the observed object (Schmidt-Kaler 



the flux at the effective filter wavelength multiplied with the 
filter bandwidth. 
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6. Discussion 

6.1. Peak Bolometric Flux 

The bolometric light curves presented in Fig. || were 
constructed from UBVRI light curves where available 
(SN 1989B, SN 1991T and SN 1994D). All others are based 
on BVRI light curves with a correction for the missing U 
band applied as described in Sect. |5.1.l| . The luminosities de- 
pend directly on the assumed distance moduli to the individ- 
ual supernovae. We have used the current best estimates from 
the literature (Table [j]), but some of the distance estimates may 
change when more accurate distances become available. We 
also corrected the magnitudes for extinction, which, for some 
objects, can be fairly substantial (SN 1989B, SN 199 IT, and 
SN 1994ae). From Table ^ it is clear that our sample displays a 
rather large range in luminosities, both in B and the bolomet- 
ric maximum. SN 1991bg is 11 times less luminous in B than 
the brightest SN la in the sample, SN 199 IT. For the bolomet- 
ric maximum we find a factor of 12.3. Excluding this peculiar 
object we still find a range in luminosity of a factor 4.7 in the 
filter passband and a factor of 3.3 in the bolometric peak flux. 
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Fig. 9. Rise to and decline from maximum in bolometric light 
curves. The time between maximum luminosity and half its 
value is plotted. The sample of SNe was extended to the one 
used in Sect. the SNe from Table |l] are shown as filled sym- 
bols. Only SNe with sufficient coverage (rise: first observation 
not later than 3 days after t_i/2> decline: first observation be- 
fore i max (Bol)) have been included. 



6.2. Light Curve Shape 

The shapes of the bolometric light curves are unaffected by the 
distance modulus and are only marginally influenced by red- 
dening (see Sect. ||). The most striking feature in these light 
curves is the inflection near 25 days past maximum (Suntzeff 



1996 ). It is observed in all SNe la, with the notable exception 
of SN 1991bg (Fig. ||). The strength of this shoulder varies 
from a rather weak flattening in the bolometric light curve of 
SN 199 IT to a very strong bump in SN 1994D. This bump 
arises from the strong secondary maximum in the R and / light 
curves. The secondary maximum is also observed in the near- 
IR JHK li ght cur ves (Elias et al. |1985| , Meikle |2000t Meikle & 
Hernandez 2000 ) which have not been used for the calculation 



of bolometric light curves presented here. 

Inspection of Fig. |] clearly shows the brighter SNe having 
wider primary peaks than the fainter SNe. 

In order to quantify this statement we have measured the 
width of the bolometric light curve at half the peak luminos- 
ity (with t_i/ 2 denoting the time it takes to rise and t+1/2 to 
decline, see Table ||). There are sufficient pre-maximum obser- 
vations available for six supernovae (SN 199 IT, SN 1992 A, 
SN 1992bc, SN 1992bo, SN 1994D, and SN 1994ae) from 
which the full width of the peaks can be reliably determined. 
In three cases the first observations were obtained well below 
the half-maximum flux level; in the case of SN 199 IT we ex- 
trapolated our bolometric light curve by about 1 .5 days and for 
SN 1992A and SN 1992bc by 2.5 days. The bright SN 1991T 
and SN 1992bc show a peak width of 26 and 23 days, re- 
spectively, the intermediate SN 1994ae one of 22 days, while 
the fainter SN 1994D, SN 1992A, and SN 1992bo remained 
brighter than half their maximum luminosity for about 18 to 19 
days. 



The pre-maximum rise t_ l / 2 in the bolometric light curve 
is, in all cases, substantially faster than the decline by 
between 20 to 30% or a time difference of between 2 to almost 
4 days (Table §). 

The decline to half the supernova's peak luminosity i+1/2 
varies from 9 days to 14 days. A weak correlation between 
maximum brightness and decline rate can be seen (Fig. [)]). The 
sample of SNe was extended to the one used in Sect. || In Fig. ^ 
only those SNe with sufficient time coverage are shown. 

Significant coverage of the rise to maximum is available for 
only two supernovae. There are too few SNe to make a defini- 
tive statement here about the rise times of SNe la in general. 
Nevertheless, our formal fits to SN 1994D and SN 1994ae give 
rise times of 16.4 and 18.2 days, respectively, for their bolo- 
metric light curves. The formal errors in this para mete r, derived 
with the Monte Carlo algorithm described in Sect. 3.2, is 2 days 
for SN 1994D and 1 day for SN 1994ae. If the bolometric lumi- 
nosities are calculated by fitting the bolometric light curve from 
the individual observations, the derived rise times are slightly 
different (18 and 19 days, respectively). A previous analysis of 
SN 1994D derived a rise time for th e bolometric light curve of 
18 days (Vacca & Leibundgut 1996 ) in fair agreement with our 
current analysis. 

These rise times are significantly larger than most light 
curve calculation of explosion models (Hoflich et al. [1996 , 



Pinto & Eastman 2000a[ ). These theoretical calculations typi- 
cally yield rise times of 15 days. 

The secondary bump clearly indicates a change in the emis- 
sion mechanism of SNe la. Its varying strength might be related 
to the release of photons "stored" i n the ej ecta (the "ol d" pho - 
tons described in Pinto & Eastman 2000b and Eastman 1997). 



In such a scenario the bump arises from the fact that the ejecta 
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Table 2. Absolute B magnitudes and bolometric luminosities. The bolometric luminosities have been corrected for the missing 
U band where appropriate. The nickel mass is derived from the luminosity for a rise time of 17 days to the bolometric peak. 



SN M B Amg logL bo i Amf 5 °' M Ni t_ 1/2 t +1/2 

(mag) (mag) (ergs -1 ) (mag) (Mq) (days) (days) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


SN1989B 


-19.37 


1.20 


43.06 


0.91 


0.57 




13.1 


SN1991T 


-20.06 


0.97 


43.36 


0.83 


1.14 


11.6 


14.0 


SN1991bg 


-16.78 


1.85 


42.32 


1.42 


0.11 




8.8 


SN1992A 


-18.80 


1.33 


42.88 


1.15 


0.39 


8.6 


10.3 


SN1992bc 


-19.72 


0.87 


43.22 


0.93 


0.84 


10.1 


13.0 


SN1992bo 


-18.89 


1.73 


42.91 


1.27 


0.41 


8.3 


9.6 


SN1994D 


-18.91 


1.46 


42.91 


1.16 


0.41 


7.5 


10.8 


SN1994ae 


-19.24 


0.95 


43.04 


0.97 


0.55 


9.6 


12.6 


SN1995D 


-19.66 


0.98 


43.19 


1.00 


0.77 




12.2 



become optically thin during this phase, due to subtle differ- 
ences in the ejecta structure and ejecta velocity. 

The exponential decline of the bolometric light curves be- 
tween 50 and 80 days past bolometric maximum is remark- 
ably similar for all supernovae. The decline at this phase is still 
dominated by the 56 Co decay and the 7— ray escape fraction 
(Milne et al. 1999). We find a decline rate of 2.6±0.2 mag per 
100 days. The sample is still very small but SN 1991bg stands 
out with a decline rate of 3.0 mag/(100 days). The uniformity 
of the late declines indicates that the differences observed at 
earlier phases are due to photospheric effects when the opti- 
cal depth for optical radiation is large and not because of the 
explosion mechanism. In particular the fraction of energy con- 
verted to optical radiation at late phases appears to change in 
an identical fashion for the majority of the objects, although 
SN 1991bg proves to be the exception to the rule once again. 
The change in the column density in the different supernovae 
must be very similar, despite the different luminosities at these 
epochs. This indicates that the kinetic energy is somehow cou- 
pled to the pre-supernova mass. 

The decline rate during the same period in the V filter is 
identical with 2.6±0.2 mag/(100 days) as is also confirmed by 
the constant fraction of the V to the bolometric flux (cf. Fig. ||). 
These decline rates have been shown to correlate very well with 
measurements derived by others (Vacca & Leibundgut 1997[ ). 
The comparison with other determinations of the decline rate 
is difficult as the slope of the decline continues to vary even at 
late epochs (e.g. Suntzeff [199^ , Turatto et al. |1996| ). This de- 
pendence on the epochs observed has to be considered when 
a comparison is attempted. In most cases the decline was es- 
timated out to phases of 200 days (Turatto et al. 1990), which 
leads to smaller decline rates than we find here. In the case 
of SN 1991bg the decline rate was measured at early phases, 
which led to a significantly steeper decline estimate (Filip- 
penko et al. |1992|). 



Note that the phase range of bolometric light curves pre- 
sented in this paper is well before any ef fects of positron escape 
could be measured (Milne et al. 1999 ) or the IR catastrophe 
takes place (Fransson et al. 1996 ). 



6.3. Nickel mass 

Given the peak luminosity it is straight forward to derive the 
nickel mass which powers the supernova emission. Near max- 
imum light the photon escape equals the instantaneous energy 
input and is directly related to t he tot al amount of 5 6 Ni s ynthe- 
sized in the explos ion (Arnett 1982 Arnett et al. 1985 Pinto 
& Eastman 2000a ). We have used our bolometric peak lumi- 
nosities to derive the nickel masses (Table ||). Since we are not 
sampling all the emerging energy from the SNe la, but are re- 
stricted to the optical fluxe s, we are underestimating the total 
luminosity. Suntzeff ( 1996| ) estimated that about 10% are not 
accounted by the optical filters near maximum light. All masses 
thus have to be increased by a factor 1.1. Another uncertainty is 
the exact rise time, which is an important parameter in the cal- 
culation. We have assumed a rise time of 17 days to the bolo- 
metric maximum for all supernovae. It is likely that there are 
significant differences in the rise times and this would alter the 
estimates for the Ni mass. A longer rise time would imply a 
larger nickel mass for a given measured luminosity. Decreas- 
ing the rise time to 12 days yields only 70% Mni of the values 
given in Table |[ Such a short rise time is excluded for most of 
the SNe la, where observations as early as 14 days have been 



recorded (SN 1990N: Leibundgut et al. 1991a, Lira et al. 1998 



SN 1994D: Vacca & Leibundgut |199J), but could still be feasi- 
ble for SN 1991bg. For a more realistic range of 16 to 20 days 
between explosion and bolometric maximum the nickel mass 
would change by only q=10% from the values provided here. 
Clearly, the dominant uncertainty in the determination of the 
nickel mass stems from the uncertainties in the distances and 
the extinction corrections. 

For a few of the supernovae, nickel masses have been mea- 
sured by other methods. SN 1991T has an upper limit for the 
radioactive nickel produced in the explosion of about 1 Mq 
based on the 1.644/im Fe lines (Spyromilio et al. 1992). This 



value depends on the exact ionization structure of the super- 
nova one year after explosion and a conservative range of 0.4 
to 1 Mq had been de rived. This is fully consistent with our es- 
timate. Bowers et al. ( 1997 ) have derived nickel masses for sev- 



eral SNe la in a similar way. Their best estimates for SN 1 99 1 T, 
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SN 1994ae, and SN 1995D are all about half the value found 
here, when converted to their distances and extinctions. How- 
ever, they point out that their values should be increased by a 
factor of 1.2 to 1.7 to account for ionization states not included 
in their analysis. With this correction we find a good agree- 
ment. Cappellaro et al. ( 1997| ) derived masses from the late V 
light curves. There are four objects in common with our study: 
SN 1991bg, SN 1991T, SN 1992A, and SN 1994D. Adjusting 
the determinations to the same distances and re-normalizing to 
our SN 199 IT Ni mass, we find a general agreement, although 
there are differences at the 0. 1 Mq level. 

Nickel masses were also derived from the line profiles of 



shapes indicates subtle variations in the energy release of these 
explosions. Pinto & Eastman (2000b) try to explain these sec- 



[Fe II] and [Fe III] lines in the optical by Mazzali et al. (1998) 



These measurements depend critically on the ionization struc- 
ture in the ejecta and had been normalized to a nickel mass of 
SN 199 IT of 1 Mq. When we scale their masses to our mea- 
surement we find a reasonable agreement. 

The nickel masses derived in our analysis are well within 
the bounds of the current models for SN la explosions (Hoflich 



etal. 1996, Woosley & Weaver 1994). There seems to be no real 
difference in 56 Ni produced by the various explosion models 
and hence a distinction by the bolometric light curve alone is 
not possible. 



7. Conclusions 

Fitting light curves with a descriptive model has a number of 
advantages over template methods. It is suited to explore the 
variety among SN la and provides an independent way to look 
for correlations. A simple application of our fitting procedure 
has been presented to demonstrate the complicated nature of 
SN la emission through the occurrence of the peak luminos- 
ity in individual filters. The most obvious signatures for non- 
thermal emission from SNe la have so far been the infrared 



light curves (Elias et al. 1985) and the lack of emission near 
1.2^im (Spyromilio et al. |1994|, Wheeler et al. |1998[). It had 



also been demonstrated through spectral synthesis calculations 



(Hoflich et al. 1996, Eastman 1997). Although there is a fairly 



large scatter of about 2 days in the relative epoch of filter max- 
imum light a clear trend to earlier maxima in R and / is ob- 
served. In fact, the / maximum occurs clearly before the B 
maximum, a trend also observed in the JHK light curves (Elias 



et al. 1985, Meikle 2000) 



Another application of the continuous approximation of 
the observations is the construction of bolometric light curves. 
Bolometric light curves form an important link between the 
explosion models and the radiation transport calculations for 
SNe la ejecta. We have demonstrated that the effects of missing 
passbands, distance modulus, reddening, light curve fitting and 
the integration methods are not critical to the shape. Amongst 
these, the uncertainties from the distance modulus towards in- 
dividual supernovae dominate. 

The shapes of the bolometric light curve of individual 
SNe la vary significantly. The secondary maxima observed in 
the R and / light curves show up with varying strength in 
the bolometric light curves as well. The variety of light curve 



ondary maxima as due to the distribution of material in the in- 
ner core of the explosion, possibly connected to the explosion 
mechanism. If this is the case, then the detailed study of the 
bolometric light curves and their differences among individual 
SNe la might provide a direct "window" into the explosion. 

With the currently best available distances we find that the 
peak luminosities of the SNe la in our sample display a rather 
large range. They imply a factor of 2.5 in the 56 Ni masses. 
SN 1991bg produced about 8 times less 56 Ni than the bright- 
est object, SN 1992bc. The range of nickel masses indicates 
significant differences in the explosions of SNe la. From the 
late-phase decline, we find that the change of the decline rate 
is rather uniform indicating similarity in change of the 7— ray 
escape fraction for all SNe la independent of the amount of 
nickel produced in the explosion. 

Acknowledgements. We are grateful to Leon Lucy for help with the 
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